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ABSTRACT 


A computer program was developed to obtain the wind 
tunnel wall corrections for wing angle of attack, induced 
drag, and pitching moment in incompressible flow. The 
vortex lattice method is used for computation of these 
correction factors. The program can be applied to wind 
tunnels of arbitrary cross-sectional shape, and wings of 
any desired planform, subject to the constraint of straight 


leading and trailing edges. 
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Riot. OF SYMBOLS 
Wing aerodynamic influence coefficient matrix 
aGemeneme, /9 Chord position 
Wing aspect ratio = b*/s 
Wing Span 
Wind tunnel cross-sectional area 
Local wing chord, streamwise direction 
Average wing chord, (s/b) 
Wang Span loading,™ (2/q) 
Wing lift coefficient = Lift/gqS, dimensionless 


Ground effect ratio 


Pitching moment coefficient = Moment/qSc 
Reet. chord 
hips. chord 


Wing induced drag coefficient 


Spanwise induced drag distribution (d/q) 
yang wnaweea drag distribution, drag/unit span 


Wing aerodynamic influence coefficient matrix 
at the wing quarter-chord 


Wall vortex loop strength 

Wing height above the ground 

Wind tunnel height 

Wing running span load, lift/unit span 
Number of horseshoe vortices per semispan 


Dynamic pressure 





Vpr, 


Suleseripts: 


FA 


Wall aerodynamic influence coefficient matrix 
Influence coefficient matrix for wall-on-tail 
Influence coefficient matrix for wall on wing 
Wing area 

Wing tape ratio, Cp/Cp 

Freestream velocity 

Induced downwash velocity 

Wind tunnel width 

Induced downwash angle at 0.75 chord position 
Influence coefficient matrix for wing-on-wall 
Inglu@ence coefficient matrix for wing-on-tail 
Axial coordinate 

Vertical coordinate 

Spanwise coordinate 

Angle of attack 

Wind tunnel wall interference factor 

Vortex strength 

Induced downwash angle at the tail 

Wing spanwise stations (5) 


foablverrieciency factor 


free air 
ground 
Wee vonreox LOGS 


wind tunnel 








TL ava. 1 
W wing 


( ) one) / 0 


Matrix Notation: 


ff Rectangular matrix 

@ } Column matrix 

: J Diagonal matrix 

[ oe Inverse of a square matrix 
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I. INTRODUCTION 


In the process of wind-tunnel testing the experimental 
procedure itself constrains and alters the flow past the 
wing under examination. This modified flow induces errors 
which must be corrected if accurate results are to be 
obtained. The major interference effects in subsonic wind 
tunnels result from the constraining influence of wind 
tunnel walls, disturbances from model supports and measuring 
equipment, and turbulence and non-uniformity in the airflow 
meselr. 

Russell [Ref. 10] treats the problems resulting from 
model support interference and introduces a method for 
reduction of this effect. The differential equations 
governing flow within a wind tunnel are the same as those 
in free air. However, the presence of the wind tunnel walls 
changes the outer boundary conditions. The result is a 
series of perturbations of the velocity-potential gradient 
normal to the model surface. The effect is pervasive in 
that almost all of the measured aerodynamic quantities must 
be corrected. 

Generally, Prandtl is acknowledged as having established 
the foundations for research on wind-tunnel wall interfer- 
ence. His lifting-line theory and concept of trailing vor- 
tices opened the doors to much experimental investigation 


in this area. The resulting studies considered two-dimensional 


eZ 








and three-dimensional effects, open and closed wind tunnels, 
and various geometrically shaped test sections. 

References 1 and 8 offer a comprehensive treatment of 
wind-tunnel wall interferences. Much tabular and graphical 
data are presented for different tunnel shapes and wing 
Spans. Pankhurst and Holder included much data on octagonal 
wind tunnels, while Pope concentrated more on a uniformly 
loaded wing. 

With the advent of sweptback, slender wings designed 
for high-speed flight, the lifting line model was no longer 
Suitable. In these situations, the wing lifting system 
may be represented by distributed horseshoe vortex elements 
(Refs. 6 and 7]. In Ref. 9, data, graphs, and formulas are 
easembled from various sources and systematically analyzed 
for various types of wings in wind tunnels. Although the 
mathematical formulations are tedious, the resulting graphs 
and tables summarized very well the data obtained from each 
theoretical treatment. 

Even with all this information available, it is still 
an arduous task to approach tabulated data with a specific 
wind tunnel shape and geometric configuration, and hope for 
a solution to the wall interference problem. It is the 
purpose of this paper to develop a computerized, closed- 
form solution to this problem. The design goal was to 
develop a digital computer code that would be applicable to 


wind tunnels of arbitrary cross-section and wings with 
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varying planforms. Such versatility would enhance the 
accuracy of test data and provide a valuable, accurate 


tool in wind-tunnel testing. 
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II. METHODS USED 


A. GENERAL 
Wind tunnels provide an invaluable tool in the design 
and performance prediction of wings and wings-and-bodies. 
However, the wind-tunnel testing environment introduces 
its own set of problems. The interference action from the 
walls of the wind tunnel alters the development of induced 
velocities by the wing lifting system. As a result, downwash 
formation iS suppressed and the wing aerodynamic character- 
istics are those of a wing having greater aspect ratio. The 
accurate determination of wind tunnel wall interference is, 
therefore, an important consideration in wing experiments. 
mae classieal wall correction factor, 6, for the case 
of a single horseshoe vortex representing the complete wing 


is defined by: 


Cc 
Q 
T 
On 
QIN 
o 
Fo 


where 
W 
Aa ~ tan V 


Using small angle approximations, the tangent is approxi- 
mately equal to the angle. The wall-correction factor may 
be expressed in terms of wing circulation, vortex span and 


average induced downwash velocity by: 


Vis) 





fee a 
WwW 


Calculation of 6 is dependent on the determination of 
the downwash velocity induced by the wind tunnel walls. 
Classical methods involve the utilization of an image 
system of vortices which cancel the effect of the wing 
lifting system at the wall. The boundary condition of no 
flow through the wall is, therefore, preserved. This system 
can Provide an accurate, closed-form solution to the ground 
effect problem. However, when four walls are present, each 
wall must be represented by an image system which affects 
other image systems, as well as the model itself. The 
result is a series of infinite systems. Truncation of this 
system introduces inaccuracies of unknown magnitude and 
Signiltilcance. 

This method of images has been applied to tunnels of 
rectangular cross-section with success. However, it is 
inadequate when arbitrary tunnel cross-sections are 
encountered. Add to this a wing with sweep and taper, and 
the problem becomes significantly altered. 

Therefore, a new method was sought to solve this problem 
Speetetedieaetor the NPS 3.5 X 5 ft. octagonal wind tunnel. 
The technique chosen was presented by Joppa in Refs. 4 and 
5. This approach lends itself very well to computer solution 
and can be applied to wind tunnels of any cross-sectional 


shape. 
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ois steway, an Extension Of JOoppa's vortex lattice 
method was made to allow the treatment of wings with varying 
Shape, sweep, and aspect ratio. Such treatment requires a 
modification of the single horseshoe vortex representation 
of the wing lifting system. The method of Weissinger as 
presented in Ref. 6 provides an accurate means of representing 


variable planforms. 


bowe vernon. LATTECE METHOD 

Since the method of images was considered to be limited 
in application, it was necessary to find another, more versa- 
tile method of representing the wind-tunnel walls. The 
problem is essentially one of finding a vorticity distribu- 
tion which will accurately model the wind-tunnel walls. MThis 
vortex system must satisfy the boundary condition of no 
flow through the walls. 

In Ref. 4, Joppa presents a method which accurately 
models the wind-tunnel walls. This was accomplished by 
representing the tunnel walls with a lattice network of 
rectangular vortex loops. This vortex system lies in the 
plane of the tunnel walls and satisfies Helmholtz's theorem 
that a vortex filament can neither begin nor end in the 
flow. Fig. 1 shows a representation of this system as 
Sooltecdmeomitne NPS 3.5 by 5 £t. wind tunnel. 

Each vortex segment in a given rectangle has vortex 


strength Daa eudmias a COrEresponding control point located 


Jey 








in the center of each rectangle where the no-flow boundary 
condition is satisfied. The last vortex lattice rectangles 
located far downstream are represented by horseshoe vortices 
trailing to infinity. The reasoning is that far downstream, 
where the influence upon the wall vortex lattice array of 

the bound vortex elements on the wing is small relative to 
the trailing vortex elements, one finds that the vortex 
strengths of the lattice terms do not change in the stream- 
wise direction. Consequently, by a straightforward cancella- 
tion argument, one can visualize replacing the vortex lattice 
rectangles by horseshoe vortices trailing to infinity. 

The mathematical presentation of the theory is based on 
the Biot-Savart law describing the velocity induced at a 
point by a vortex segment. Joppa in Refs. 4 and 5 follows 
the mathematical formulation beginning WilEn the Biot=<SavaEt 
law. This development is summarized in Appendix B. 

The general matrix representation of this method centers 
Aeeuncecne boundary condition of no flow through the wall. 
As mentioned previously, this condition is satisfied ata 
control point located at the center of each vortex loop. 
Maeretore, it follows that the velocity induced at each 
control point by the vortex loops must be exactly cancelled 
by the velocity created by the wing lifting system. This 


is presented below in matrix representation 


[ROL] i = Von! ers (3) 


Ie, 





which allows a solution of vortex lattice strengths in 
terms of the wing vortex strengths. 


(10) [ROL] 


[WoL] iT} (4) 
Meee that the [ROL)] matrix, although possibly quite large, 
1s square and hence may be inverted. 

[ROL] represents the wall influence coefficient matrix. 
It is obtained by calculating the downwash induced at each 
wall control point by a unit value of each vortex loop. 
Therefore, Beles 1s the downwash at wall control point "i" 
Eioawecd by vortex Loop "j". [WOL] is the influence coeffi- 
client matrix describing the wing effect on each loop con- 
tao Doint. Blas is the downwash produced at control point 
[= oy the ae wing horseshoe vortex. 

By inverting [ROL] and premultiplying with [WOL], the 
strength of the wall vortex loops can be found in terms 
of the wing circulation. The inversion of [ROL] can be 
very time-consuming, since it may bea very large matrix. 
However, by using symmetry, the order of this matrix can 
be reduced by a factor of four. This requires the wing to 
be placed on the horizontal and vertical centerline of the 


wind tunnel. 


Cee Reo REPRESENTATION 
iieemn Ss, Opa Stated that the wing can be repre- 


sented by a lifting system more complicated than a single 
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elementary horseshoe vortex. In order to model swept wings 
of arbitrary planform, a more definitive lifting system 
model was employed. Gray and Schenk, Ref. 6, present a 
modified Weissinger method designed to find the load dis- 
tribution of wings with arbitrary planform. The wing is 
represented by a series of horseshoe vortices. The midspan 
of each horseshoe vortex bound element is located at the 
heucmauarter—-chord position. The condition of flow tangency 
is satisfied at the wing three-quarter-chord point in the 
center of each horseshoe vortex. Figure 2 illustrates the 
Senet gGuracron chosen. 

Pearson [{Ref. 7] modified this technique by allowing 
the bound vortices to be swept, and including chord-wise 
horseshoe vortices. For this presentation it was felt that 
unswept bound vortices would accurately represent the lifting 
system of a swept wing. In addition it was felt that little 
deterioration of accuracy would result if the chordwise 
system of bound vortices was combined into a single bound 
Mercex at Eme quarter—chord. Figure 2 illustrates the 
@ontiguraeron chosen. 

An important concept in this representation of the wing 
is that of flow tangency at the three quarter-chord. The 
downwash angle at the wing can, therefore, be determined 
Eyota maga concrol point at this position in the center 
of each horseshoe vortex. The rationale behind the three- 


quarter chord flow tangency condition is a product of 


Zk 





Figure 2. 


Wing Geometry 


OAPs 








two-dimensional flow considerations. Extension of this 
technique to the three-dimensional wing has been shown to 
produce experimentally confirmed results. However, it is 
still unclear why this technique works so well. 

The wing aerodynamic influence coefficient matrix was 
then developed. [AERO] represents the downwash induced at 
each 0.75 chord wing control point by a unit value of wing 
horseshoe vortices. The matrix formulation for induced 


downwash formation in free air 1S presented below: 
{w} = [AERO] iT, } (5) 


Dividing the above expression by freestream velocity, one 
obtains an expression for 0.75 chord control point induced 
downwash angles. The flow tangency requirement equates 
this column matrix to the wing angle of attack column 


vector. 
W = r 
{=} = [AERO] Sa (6) 


and 


30/4 


. Zh 
It should be noted here that = 1s equal to TT from 
the Kutta-Joukowski Law. Substitution into equation 6 


yields the following results: 


ZS 








(¥} = 0.5 [AERO] {2} (7) 
Vv dpa 


Solving equation 7 for (=}, one can obtain the wing span 
Head aiStribution. Since = is equivalent to the sectional 
value of Cyc, integration of the = distribution vs. span 


provides a direct measure of the wing lift coefficient, C 


L° 
When one assumes a uniform angle of attack value of one 
radian along the span, the resultant integration of the 
Span load solution will yield the wing lift curve slope, 
Cc... It should be borne in mind that these techniques are 


La 


predicated upon principles from thin airfoil theory and 
in the limit for an infinite aspect ratio wing will yield 


a wing lift curve slope of 2u(raa +). 


De DETERMINATION OF THE WALL CORRECTIONS 

The classical wall interference factor 6, which includes 
a ratio of wing to tunnel cross-section area, is generally 
found by representing the wing as a single horseshoe vortex. 
Since the present formulation utilizes ten horseshoe vor- 
tices per semi-span, § in a modification can best be found 


by continuing our matrix algebra. 





on = Cn + oC, 
tna = @ “L (8) 
a 
FA 
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Syeocsumings that the lift. coefficient is the same in 


both tunnel and free air, the following equations result: 








+ 6 Cc. (9) 








(10) 


The free air lift slope was found in the previous section. 
The wing lift curve slope for the wing in the wind tunnel 
can be found by a simple extension of this same method. 

The total downwash at the wing control points is equal to 
that induced by the wing added to that created by the tunnel 
walls. The effect produced by the tunnel walls is repre- 
sented by the vortex rectangle-on-wing influence coefficient 
matrix [ROW]. re is the downwash induced at the yeh wing 
Cemtrol point by the a wall vortex ooo. Multiplying this 
matrix by the vortex loop strength, one obtains the total 


wall effect. 


{w} = [AERO] oe + [ROW] (eo ia) 


Ps: 





Substitution of equation 4 for 1 forms the following 
set of equations: 


iL 


{w} = [AERO] {tT} + [ROW] [ROL] ~ [WOL] ee 


} = 0.5 [{AERO] + [ROW] [ROL] ~+ {wor} {4} (12) 
3 . 
4 

By establishing one radian angle of attack, solving 
fOr = and integrating the result, the wind-tunnel lift 
curve slope can be determined. It should be noted that the 
correction factor obtained from this process is different 
from the classical 46, since the correction factor computed 
above already contains the wing area to tunnel area ratio 
influence. Therefore, the correction to angle of attack 


can be obtained by a direct application of this 64. 
cor — ero C (E33) 


In most literature, the correction for induced drag 
is obtained by application of the same 6 as that used for 


angle of attack. 


AC Se ae (14) 
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However, in this presentation, the induced drag correc- 
tion is calculated in a different manner because of added 
available aerodynamic information. The spanwise induced 
drag distribution in free air is given by the following 


Formula : 


f=) [pacRO] {—} (15) 
eS Tap a FA q _ 


[DAERO] is the wing aerodynamic influence coefficient 
matrix computed with the wing control point now located 
at the gquarter-chord position. The wall effect can similarly 


be found by including additional matrix terms for the wall 


effect on the wing quarter-chord position. 


{=} - Pew [[DAERO] + [ROW] [ROL] ~~ (WOL] ] {=} (16) 


T As Au 


imtegmarion Ofsequations 15 andwlGé resuit in,the free 
air and wind tunnel induced drag, respectively. The induced 
digeaq cormection factor can then be found in the following 


manner. 


2 
5 7 REA ey 
tPa 
2 
“Da ie OE ae 
T 


a 





pawn op. Prk. 
tra 
: 2 
Kn ss oF [Cy 
1p lb 
(Eze) 
3 = Sey Sn 
2 
Ce 


When testing a model in the tail-on configuration, 
the presence of the wind tunnel walis serves to modify the 
downwash angle induced at the tail. This altered flow 
field can have a significant effect on the pitching moment 
coefficient. The following increment in pitching moment 


coefficient is generated by the horizontal tail: 


ae ee Sa. ya (18) 


The term in parentheses 1S approximately equal to the 
21 
Stabilizer effectiveness, x5 Weevil teoandgleoof attack 
at 


can be estimated as follows: 
o} =-C@=—- € + 6. (19) 
ak 


At a given value of wing lift coefficient an incremental 
change in downwash angle at the tail due to the wall effect 


can be computed: 
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(20) 


Giererore, the correction to pitching moment coefficient is 


given by the following relation: 


(21) 


Since the downwash angle € is equal to * at the tail 


position, « for free air and for the wind tunnel can be 


hound trom the £Ollowing relations, respectively: 


3 & 
5 = Q.5[WOT] ee 


FA FA 


€ = (223) 


7 O.S5[[WoT] + [ROT] [GAMAL] ] (3 
qt 


[WOT] is the wing effect on the tail position and 
[ROT] [GAMAL] represents the wall effect on the tail. 


Equation 20 can then be used to find Ae. The pitching 


moment correction, as evaluated in equation 21, is thus 


dependent on being given an experimental value for the 


tail effectiveness parameter. 
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tout oe AND! DESCUSSION 


A. COMPUTER PROGRAM 

The computer program developed consists of two major 
Subroutines. Subroutine Loop calculates the downwash at 
any point due to a vortex rectangle. Subroutine HSHOE 
computes the downwash at any arbitrary point due toa 
horseshoe vortex. HSHOE is used for both the wing horseshoe 
and those of the wall system located far downstream in the 
tunnel walls. The program offers two basic options. One 
option establishes whether ground effect or wind-tunnel 
wall corrections are to be calculated. The other option 
indicates whether the tail-on or tail-off case is to be 
considered. 

Each computer run offers the inclusion of many wings 
of different span, aspect ratio, taper ratio, and sweep. 
The wing parameters are read as input data after the inver- 
sion of the large wall-on-wall coefficient matrix. There- 
fore, with a given tunnel configuration, a variety of 
different wings can be included with a minimum of computer 
time consumed. 

Integration of the wing span loads in accomplished by 
a subroutine obtained from the NPS source library. It 
should be noted that in the integration scheme the assumption 
is made that the wing loading at n = 0.0 is the same as that 


for n = .05b. Any inaccuracy resulting from this assumption 
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appears to be negligible. Appendix C describes the necessary 


Page tO the program. 


Be GROUND EFFECT 

Originally, ground effect calculations were conducted 
merely to check the accuracy of the vortex lattice method. 
As previously stated, the method of images provides a 
closed-form solution to the ground effect problem. There- 
fore, this provided an excellent opportunity for a conver- 
gence check. As an object of academic interest an option 
for ground effect computation was retained in the final 
CGemputer program. 

The results of the convergence check showed that the 
vortex lattice method could very accurately compute the 
velocity induced at’ a point by the presence of a wall. By 
applying the theory previously developed, the lift curve 
Slopes for free air and ground effect were computed. Results 
were then compared with those presented in Etkin [Ref. 3]. 
The free air lift curve slope was found to correlate extremely 
well for both straight and sweptback wings. 

Furthermore, in Ref. 3, Etkin defines a ground effect 
ratio as the ground effect lift curve slope divided by that 
in free air. The results are tabulated in Table l, and 
Figure 3 compares the results obtained from the computer 
program with those presented in Etkin. As can be seen, 


agreement 1s very favorable. 


con 








GROUND EFFECT RATIO 
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1.30 AR 8,0 
TR 1,0 
Sweep 

110 A Span 
AR 6.0 
TR 0.3 
Sweep 

1.09 

1.07 

}.05 

1.03 


Figure 3. Ground Effect Ratio Plotted as a Function 


of Wing Height Above the Ground 
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Several computer runs were conducted at different wing 
heights above the ground. For each run, sweep was varied 
from 0 to 40 degrees. Results are tabulated in Figure 4. 
meeean be noted that the ground effect ratio initially 
increases with sweepback but returns to the unswept value 
at about 30 degrees. 

As a note of caution, care should be exercised when 
the wing gets very close to the ground. It is necessary 
to greatly decrease the vortex loop size for accurate 
results. This procedure seems to provide a much more finely 
defined ground representation. Otherwise, calculations 


Will result in inordinately high values. 


C. WIND TUNNEL APPLICATION 

A convergence check of the wind tunnel wall correction 
factor was obtained by applying the program to a 3 X 5 ft 
rectangular tunnel. The correction factor obtained was 
converted to the classical 6 by including the factors z 
and 57.3 (radians to degrees). The result of 0.1216 was 
compared to that presented by Joppa in Ref. 5. The wing 
geometric configuration presented in Ref. 5 consisted of a 
Single horseshoe vortex. The value of 0.127 from Ref. 5 
correlates very well with the results of this study. 

Tables 2 and 3 present the data obtained from compu- 


eaetons with the NPS 3.5 X 5.0 £t octagonal wind tunnel. 


Computer runs were conducted using wings of varying span 
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Figure 4. Wing Span Load Distribution 
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and sweep angle. Figures 4 and 5 illustrate the span loads 
eemedtared for a typical wing. It is interesting to note 
that as sweep angle increased, negative induced drag appears 
on the outboard wing stations. Therefore, total induced 
drag is seen to steadily decrease with increasing sweep 
angle. 

As sweep angle is increased, the free air lift curve 
Slope decreases. Correspondingly, the corrections to 
angle of attack and induced drag are reduced. Figure 6 
graphically displays the data for a straight wing. In the 
Paste the Correction factors for the NPS 3.5 X 5.0 ft 
octagonal wind tunnel were obtained by assuming an elliptic 
cross-section. The correction factors presented in Ref. 1 
Her Clliptic cross-section are 0.610 for angle of attack 
and 0.0106 for induced drag. In comparison, the vortex 
lattice program computes values of 0.7086 and 0.01075, 
mesmectively. It is interesting to note that the induced 
drag correction is almost identical, whereas the correction 
Eeomangle or attack 1s slightly higher for the octagonal 
wind tunnel. Figure 7 illustrates the angle of attack 
eeorrection factor with varying tunnel span ratios. The 
correction factor can be seen to steadily increase with the 
SOan ratio. 

It is an interesting observation that the value of 
induced drag for a wing span of 5.0 ft was almost zero. 


This compares very well with airfoil theory in that a wing 
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which spans the wind tunnel represents an approximation to 
the two-dimensional wing. Induced drag for a two-dimensional 
wing should equal zero. In addition, the wind-tunnel lift 
curve slope is within 0.0006 of the theoretical two-dimensional 
Value of 2T rad > 

A sample computer run was made for the moment correction 
factor. The tail quarter-chord was placed two feet behind 
the wing bound vortex, on the tunnel centerline. The 
correction factor obtained for a straight wing of aspect 
ratio 6.0 was -0.2297. For a wing with sweep angle of 
30 degrees and taper of 0.3, a value of -0.1887 was calcu- 
lated. The above correction factors must be multiplied 
by the stabilizer effectiveness parameter, using consistent 
angular measures, in order to find the pitching moment 


correction. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


It is felt that the vortex lattice computer program 
as presented in this paper offers a very reliable and 
flexible tool for the determination of wind-tunnel wall 
GGercctions. It is capable of treating wings of any span, 
taper ratio, aspect ratio, and quarter-chord sweep, subject 
to the constraint of straight leading and trailing edges. 
The geometric shape of the wind tunnel cross-section can 
Similarly be varied at the discretion of the program user. 
Resultant accuracy appears to be excellent. However, cau- 
tion should be exercised in determining the vortex loop size. 
As the wing is positioned closer to the walls, a smaller, 
more definitive vortex loop representation must be used. 

A logical extension of this program would be the inclu- 
sion of wings with an antisymmetric load distribution. 
Furthermore, slight modifications could incorporate the 
ability to consider the wall effects on other stability 
derivatives. In its present form, this program offers 
the wind tunnel user a versatile tool for determination of 


wima tunmel wall corrections. 
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APPENDIX A 


The general assumptions used in the compilation of 


this program are listed below. 


1. The two-dimensional lift-curve slope is 27 rad ?. 


2. The wing load distribution is symmetrical. 

3. Aeroelastic effects were not considered. 

4. The wing is located on the horizontal and vertical 
centerline of the wind tunnel. 

5. The leading and trailing edges of the wing are 


Sticke iit. 
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APPENDIX B 


DEVELOPMENT OF INDUCED VELOCITY EQUATIONS 


The equations required to solve for induced velocities 
consist of two main types. One equation solves for the 
velocity induced by an arbitrarily oriented vortex segment, 
whereas the other basic equation applies to horseshoe vortex 
elements. These equations have been incorporated in sub- 
routines LOOP and HSHOE, respectively. The development of 
these equations presented below follows closely that found 
mueeret. 4 and 5. 

According to the Biot-Savart law, the velocity induced 
pea POINt p 1S found from: 
[ 


—— (cos B 


es + cos B.) V (1B) 


= | 
I 


ll 


For vortex segments oriented in the Y-Z plane, 


ee ae IR 
S24 a) oe Sn. 
cos B. + cos B = TRORIS (S (Ry Ry) ) (2B) 


i iz 1% 


It is now necessary to consider v, the unit vector 
primeneestablishes direction. From Fig. 8, 1t can be seen 


femal tc 


7 (3B) 


45 








| 
| 
| 
(XPYPZP) | 
| 
{ 


Figure 8. Velocity Induced By A Vortex Segment 
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Newamg from Fig. Sethat 


ie = |Ry X S| (4B) 
then 
i eo meee re (Rk. SS Rs) jc 
iV Z ay eZee ne ene Z 
~ TR 
e (Ry 5. Ri yS,) ] 


Calculating the induced velocity in the direction of 


outward normal of point P per unit I, 

ae ee 

rT = Pf n (5B) 
where the unit normal vector is defined by: 


oi + (sin 9p)j + (cos Op)k 


S| 
I 


Therefore, one finally obtains: 


ee ee eee ee 
87R,R, (hs) 
(6B) 
[(sin 3) OE Roe) oe SSNS) SOS ee cece | 
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A conditional statement must be included in the program 
to set the downwash equal to zero if hs is equal to zero. 
This will preclude having a zero term in the denominator. 

The trailing elements of a vortex loop can be treated 


with more Simplicity. For these segments, 


(R,? _ n2) 1/2 
cos 8 
a R, 
(7B) 
(Ree n2) 1/42 
cos B = 
2 R, 
However, it can be seen from Fig. 8 that 
: a “gp = NS + (yA - yp)? 2 (2A - 2p)? 
(8B) 
n2 = (yA - yp) + (zA- 2p) - 
Combining equations 7 and 8, 
i R 
1 
(3B) 
cos B en eee 
2 R 


As a result, the equation for the trailing elements 


of a vortex loop is simply 
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. — o—_7 misin Op) (zA = zp) 


(10B) 
+ (cos 8p) (yp - yA)] 

For the horseshoe vortex, the equation for the bound 
segment 1S identical to equation 6. Similarly, the 
trailing elements can be treated as in equation 11 with a 
Slight modification. In this case, the cosine of one 
enclosed angle is equal to 1, since the trailing element 
extends downstream to infinity. The induced velocity 


equation becomes 


1 + cos Bs 


SS 2A - 2p) 5 + (yp - ADE] 
4th 

W ie cos B 

a [ (Sine) (2auezp) 


hone 


(11B) 


+ (cos 6p) (yp - yA) ] 
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Card Number 


iL 


Input cards 5 and 6 can be repeated for different 


planforms and tail locations. 


APPENDIX C 


COMPUTER PROGRAM INPUT 


Column Number 


Description 


PrEequam t2t lve 


MM=Number of Spanwise 


Stations 


NN=Number of Downstream 


Stations 


Deltax=Loop length 


NOPT=1 for Ground 
Effect 


NTAIL=1 for tail-on 


Tunnel y-coordinates 
(one quadrant only) 


Tunnel z-coordinates 
(one quadrant only) 


span 
Aspect Ratio 
Taper Ratio 


c/4 Sweep (degrees) 
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